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Abstract. We perform a detailed forecast on how well a EuCLiD-like survey will be able 
to constrain dark energy and neutrino parameters from a combination of its cosmic shear 
power spectrum, galaxy power spectrum, and cluster mass function measurements. We find 
that the combination of these three probes vastly improves the survey's potential to measure 
the time evolution of dark energy. In terms of a dark energy figure-of-merit defined as 
(a(wo)a(w a ))~ 1 , we find a value of 454 for EuCLiD-like data combined with PLANCK-like 
measurements of the cosmic microwave background (CMB) anisotropies in a fiducial ACDM 
cosmology, a number that is quite conservative compared with existing estimates because 
of our choice of model parameter space and analysis method, but still represents a factor 
of 3 to 8 improvement over using either CMB+galaxy clustering +cosmic shear data, or 
^vj CMB+cluster mass function alone. We consider also the survey's potential to measure dark 

,/ energy perturbations in models wherein the dark energy is parameterised as a fluid with a 

nonstandard non-adiabatic sound speed, and find that in an optimistic scenario in which wq 
deviates by as much as is currently observationally allowed from —1, models with c^ = 10 -6 
and Cg = 1 can be distinguished at more than 2a significance. Under the same optimistic 
assumptions, if the Jeans mass associated with dark energy clustering falls within the cluster 
S^ mass range observed by the survey, then the order of magnitude of the dark energy sound 

speed can potentially be pinned down. Finally, we find that the sum of neutrino masses 
can be measured with a la precision of 0.01 eV, even in complex cosmological models in 
which the dark energy equation of state varies with time. The la sensitivity to the effective 
number of relativistic species N^ is approximately 0.02, meaning that the small deviation 
of 0.046 from 3 in the standard value of N^ due to non-instantaneous decoupling and finite 
temperature effects can be probed with 2a precision for the first time. 

1 Present address: School of Physics, The University of New South Wales, Sydney NSW 2052, Australia 
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1 Introduction 

The coming decade will see spectacular advances in the measurement of the large-scale struc- 
ture distribution in the universe. Perhaps the most interesting of these measurements are 
the large-scale photometric surveys to be conducted by the Large Synoptic Survey Tele- 
scope (LSST) [1] and the ESA Euclid mission [2]. Both projects will map the positions 
and measure the shapes of order a billion galaxies in a significant fraction of the current 
Hubble volume. This will in turn allow for a precision measurement of both the galaxy 
clustering and the cosmic shear power spectra, and likewise an impressively precise deter- 
mination of the cosmological parameter values. As an example, assuming a vanilla ACDM 
model extended with nonzero neutrino masses, a EuCLiD-like survey will be able to measure 
the neutrino mass sum ^ m v most optimistically at the 0.01 eV level of precision [3] when 



combined with measurements of the cosmic microwave background (CMB) anisotropics from 
the Planck mission [4]. Such a precision will see the absolute neutrino mass scale detected 
at 5ct+ confidence even if the true value of Y2 m v should be the minimum compatible with 
current neutrino oscillation data, i.e., ^ m v ~ 0.06 eV [3]. 

As yet unexplored in reference [3] is the role played by the cluster mass function in 
cosmological parameter inference. Weak gravitational lensing measurements available to 
both the LSST and Euclid will allow for the efficient detection and mass determination 
of galaxy clusters; Euclid, for example, is expected to detect and accurately measure the 
masses of close to 100 000 clusters [2]. In this work, we continue the cosmological parameter 
sensitivity forecast begun in [3] by adding the cluster mass function inferred from a Euclid- 
like cluster survey to the galaxy and the shear power spectra measurements already considered 
in reference [3]. We also extend the study to dark energy models with a time-dependent 
equation of state and/or a nonstandard non-adiabatic sound speed. As in [3] we shall adopt 
the survey specifications of the Euclid mission in terms of the number of objects observed, 
the redshift range, and the sky coverage. However, the analysis procedure can be easily 
adapted to other similar redshift surveys such as the LSST. 

The paper is structured as follows. We discuss first our dark energy parametrisation 
in section 2, before introducing in section 3 the cluster mass function as a cosmological 
observable. In section 4 we examine some uncertainties likely to be encountered in a Euclid- 
like measurement of the cluster mass function, and discuss how we model and propagate these 
uncertainties in our forecast analysis. Sections 5 and 6 outline respectively our mock data 
generation and forecast procedures, while section 7 contains our results. We conclude in 
section 8. 

2 Dark energy parameterisation 

Although dark energy is the most popular explanation for the apparent accelerated expansion 
of the universe, there is as yet no consensus on its actual physical properties. For this reason, 
and for reasons of simplicity, dark energy is usually described as a fluid obeying the laws of 
general relativity. The homogeneous part of this fluid is responsible for driving the expansion 
of the universe, and can be represented by an equation of state w{r) = Pq(t)/pq(t), where 
Pq(t) and pq(t) denote the unperturbed dark energy pressure and energy density respec- 
tively, and t is conformal time. Except in the case of a cosmological constant, for which 
w(t) is precisely the constant —1, dynamical dark energy models have in general equations 
of state that are functions of time. For this reason, we model dark energy equation of state 
using the popular parametrisation [5, 6] 

w(t) = w + w a [l - a(r)], (2.1) 

where wq and w a are constants, and a{r) denotes the scale factor. Note that this parameter- 
isation should be regarded simply as a toy model that facilitates comparisons across different 
observational probes. We make no pretence here that it actually captures the behaviour of 
any realistic dynamic dark energy model. For an example of a forecast tailored specifically 
to scalar- field models of dark energy, see, e.g., [7]. 

A general relativistic fluid evolving in an inhomogeneous spacetime will in general de- 
velop inhomogeneities of its own. We express inhomogeneities in the dark energy density in 
terms of a density contrast 5q(t,x) satisfying pq(t,x.) = pq(t)[1 + 5q(t, x)], where pq(t, x) 
is the fully time- and space-dependent dark energy density. 



Evolution of the density contrast 5q (t, x) can be described by a set of (nonlinear) fluid 
equations coupled to the (nonlinear) Einstein equation. For the nonlinear aspects of the 
formation of clusters we refer to section 3.1 and appendix A. For the purpose of calculating 
the linear matter power spectrum we implement the linear evolution for the dark energy 
density contrast, described in the synchronous gauge and in Fourier space by the equations 
of motion (see, e.g., [8, 10-12] 

5 Q + (l + w )(e Q + ^)+3(c 2 s -w)H8 Q + 9(l + w)(c 2 s -cl)H 2e ^=0, 

~2 1,2 

Q + (l- 3c 2 s )H6 Q - ^-<5 Q + k 2 a Q = 0, (2.2) 

where 5q(t, fc) now denotes the dark energy density contrast in Fourier &:-space, 0Q(r,k) 
is the divergence the dark energy velocity field, h the metric perturbation, % = a/a the 
conformal Hubble parameter, and c 2 = SPQ/5pQ\ Tes t and c 2 = Pq/pq are the non-adiabatic 
and adiabatic dark energy sound speeds respectively. Note that the non-adiabatic sound 
speed c 2 is defined as the ratio of the pressure perturbation 5Pq to the energy density 
perturbation 5pQ in the rest-frame of the dark energy fluid, while the adiabatic sound speed 
c 2 is related to the homogeneous fluid equation of state via w = 3(1 + w)(w — c 2 )l-L. 

3 The cluster mass function as a Euclid observable 

Cluster surveys can be an excellent probe of dynamical dark energy because the abundance 
of the most massive gravitationally bound objects at any one time depends strongly on both 
the growth function of the matter perturbations and the late-time expansion history of the 
universe (see, e.g., [13-18]). The Euclid mission will identify clusters in the photometric red- 
shift survey accompanied by a spectroscopic follow-up. The same survey will also determine 
the masses of the detected clusters by way of weak gravitational lensing (cosmic shear). 

3.1 Cluster mass function from theory 

A simple quantification of the cluster distribution is the cluster mass function. Denoted 
dn/dM(M, z), the cluster mass function counts the number of clusters per comoving volume 
in a given mass interval [M, M + dM] as a function of redshift z. 

For any given cosmological model, an accurate prediction of the corresponding cluster 
mass function necessitates the use of A-body/hydrodynamics simulations. However, a num- 
ber of fitting functions, calibrated against simulation results in the vanilla ACDM model 
framework, have been proposed in the literature (e.g., [19-21]). In this work, we model the 
cluster mass function after the Sheth-Tormen fitting function [19] 
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where the fitting parameters are a = 0.707, A = 0.322184, and p = 0.3, and p m (z) is the 
mean matter density (it was shown in [45] that this yields a very good fit also in models 
with non-zero neutrino mass). The quantity a^(M, z) denotes the variance of the linear 
matter density field smoothed on a comoving length scale X sm = a _1 [3Af/(4-7r / o m )] 1 ' 3 , and is 



computed from the linear matter power spectrum P^(k, z) via 

1 f'°° 
al(M,z) = ^ J dk k 2 \W(kX sm )\ 2 P^(k,z), (3.2) 

where W(x) = 3 (sinx — xcosx)/x 3 is the Fourier transform of the spherical (spatial) top-hat 
filter function. The linear power spectrum P^(k, z) can be obtained from a Boltzmann code 
such as Camb [22]. 

The quantity 6 C (M, z) is known as the linear threshold density of matter at the time 
of collapse. Its value is established by tracking the full nonlinear collapse of a spherical 
top-hat over-density, noting the time r co ii the region collapses to an infinitely dense point, 
and then computing from linear perturbation theory the linear density contrast at r = 
Teoii. In many applications it suffices to take the constant value 5 C = 1.68. In dark energy 
cosmologies, however, this may not be a very good approximation (see, e.g., [9, 23]). Here, 
we estimate 5 c (M,z) as described immediately above, and track the spherical collapse of a 
top-hat overdensity by solving the equations 

A +n^ = - 4 ^a 2 [p m 5 m + p Q (l + 3c 2 s )5 Q ], (3.3) 



XX 3 

3 



5 m (r) = [1 + S m (n)] 



X{Ti) 

X(r) 



(3.4) 



where X(t) is the comoving radius of the top-hat, <5 m (r) and 5q(t) the matter and the 
dark energy density contrasts respectively in the top-hat region, and Tj is a reference initial 
time. Note that equation (3.4) follows from conservation of the total mass of nonrelativistic 
matter M m in the top-hat region. For more detailed discussions of the spherical collapse 
model, we direct the reader to references [24-26]. 

The presence of the pn(l + 3c 2 )Sq term on the right hand side of equation (3.3) in- 
dicates that the dark energy component also participates in the collapse, especially when 
the initial dimension X(n) of the top-hat matter overdensity exceeds the comoving Jeans 
length associated with the fluid's non-adiabatic sound speed [9, 27-29]. The resulting linear 
threshold density 5 C therefore exhibits generically a dependence on the mass of the collapsing 
region, in addition to the usual z-dependence. However, tracking the nonlinear evolution of 
the dark energy density contrast 5q is in general nontrivial because the spherical top-hat 
region is well-defined strictly only in the c 2 = and the c 2 — > oo limits, where, supplemented 
with [9, 28] 

Pq(1 + 3c s 2 ),5q^0, c s 2 ^oo, (3.5) 

PQ + z(n + ^)( PQ + P Q ) = 0, c s 2 = 0, (3.6) 

the collapse equation (3.3) can be solved exactly. Extending the application of equation (3.3) 
to the intermediate regime necessitates additional assumptions, which do not however render 
the system any less intractable [9, 29]. For this reason, we shall resort to modelling the mass- 
dependence of 5 c (M,z) by interpolating between the two known limits using a hyperbolic 
tangent function, where the location of the kink at each redshift is adjusted to reflect the 
Jeans mass corresponding to the given c 2 . See appendix A for details. 

Finally, the virial radius R Yir and the virial mass M Y - ir can likewise be computed from 
equations (3.3) to (3.6) in the two limits of c 2 . Here, i? v ir = o.X(T v i T ) is defined as the physical 



radius of the top-hat region and M v j r the total mass contained therein at virialisation, where 
virialisation is taken to mean the moment at which the virial theorem is satisfied by the 
collapsing region, r v ; r . The virial mass M wir = M m + MQ(r v i r ) counts both contributions from 
nonrelativistic matter M m and from the clustered dark energy Mq = (47t/3)/9q<5q(<2X) 3 , and 
it is M v i r , not M m , that we identify with the cluster mass M throughout this work. Further 
details can be found in appendix A. The virial radius will be used in section 5.1 to determine 
the cluster mass detection threshold. 

3.2 The observable 

The actual observable quantity is the number of clusters iVj au\ in the redshift bin i and the 
(redshift-dependent) mass bin j(i), defined as 

r°° r°° A 2 v dn 

N i>m =An dz dM—-(z)W t , m (M,z) — (M,z), (3.7) 
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where Afi is the solid angle covered by the survey, d 2 V/(dildz)(z) the comoving volume 
element at redshift z, and Wjj(j)(M, z) is the window function defining the redshift and mass 
bin. Note that the window functions are in general not sharp in z- and M-space because 
of uncertainties in the redshift and the mass determinations (see sections 4.1 and 4.2). We 
defer the discussion of our binning scheme to section 5.2. 

4 Measurement errors 

4.1 Redshift uncertainty 

The Euclid photometric survey will measure redshifts with an estimated scatter ofa(z) ~ 0.03(1 + z) 
and almost no bias [2]. Nonetheless, because the detected clusters will be subject to a follow- 
up spectroscopic study, the effective uncertainty in the redshift determination per se can be 
taken as negligible. Additional redshift errors may arise from the peculiar velocities of the 
clusters, where velocities up to ~ 1000 km s _1 may lead to an error of 5z ~ 0.003. However, as 
we shall see in section 5.2, even the narrowest redshift bins adopted in our analysis typically 
have widths of order Az ~ 0.03, i.e., a factor of ten larger than the peculiar velocity uncer- 
tainty. We therefore treat the cluster redshift as infinitely well-determined, and approximate 
the window function of the redshift bin i as 

Vi(z) = 6{z - ZmiivM^max,* ~ z), (4.1) 

where z m \ a ^i and z msx ,i denote, respectively, the lower and upper boundaries of the bin. 

4.2 Uncertainty in the weak lensing mass determination 

The mass of a cluster determined through weak lensing, M b s , is subject to scatter and bias 
with respect to the true mass of the cluster M [30, 31]. For a mass determination algorithm 
that treats clusters as spherical objects, the triaxiality of realistic cluster density profiles, 
for example, could cause the cluster mass to be over- or underestimated depending on the 
orientation of the major axis in relation to the line-of-sight. Additional biases are incurred 
if the true density profile deviates from the assumed one. 



In this work, we assume that the bias can be controlled to the required level of accuracy, 
and model only the scatter in the mass determination using a log-normal distribution [30, 31], 
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Here, P (M t, s |M) denotes the probability that a cluster with true mass M is determined to 
have a mass M Q ^ S by the survey, and /x and a are, respectively, the median and the standard 
deviation of the distribution. Since we assume an unbiased mass determination, the mean of 
the distribution is M and thus fj, = In M — a 2 /2. We use a = 0.2 to model the mass scatter. 
The distribution (4.2) can be integrated over M b s in the interval [M min j(j),M max j(j)] 
in order to determine the probability that a cluster of true mass M in the redshift bin i will 
be determined to lie in the mass bin j(i). Combing the resulting integral with the redshift 
window function Vi(z) from equation (4.1), we obtain the window function for the redshift 
and mass bin {i,j(i)}, 
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where M tllT (z) is the mass detection threshold, to be discussed in section 5.1. 
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5 Mock data generation 

The observable quantity in a cluster survey is the number of clusters Niju\ in the redshift bin i 
and mass bin j(i). Thus for any given fiducial cosmology and survey specifications, one may 
compute the fiducial cluster numbers Nffa) as per equation (3.7), and then create a mock data 

set Niju\ by assuming iVj ,v^ to be a stochastic variable that follows a Poisson distribution 
with parameter Nr^,~. An ensemble of realisations may be generated by repeating the 
procedure multiple times, and parameter inference performed on each mock realisation in 
order to assess the performance of a survey. 

This is clearly a very lengthy process. However, as shown in reference [32], for the sole 
purpose of establishing a survey's sensitivity to cosmological parameters, it suffices to use 
only one mock data set in which the data points are set to be equal to the predictions of the 
fiducial model, i.e., -/Vjj(i) = Nf?,*. This much simplified procedure correctly reproduces the 
survey sensitivities in the limit of infinitely many random realisations, and is the procedure 
we adopt in our analysis. 

In the following we describe in some detail the survey specifications that go into the 
computation of the fiducial N^ (i y the mass detection threshold, our redshift and mass 
binning scheme, and the survey completeness and efficiency. In section 5.4 we summarise the 
mock data sets to be used in our parameter sensitivity forecast. 

5.1 Mass detection threshold 

We model the redshift-dependent mass detection threshold M t hr(z) following the approach 
of references [13, 33]. A cluster of mass M at redshift z produces a shear signal kg(M,z), 

where 

M/[7tR 2 (M,z)} 



kg{M,z) = a(M,z)'- 



T, CT (z) 



(5.1) 



Here, assuming a truncated Navarro- Frenk- White (NFW) density profile [34], R s (M,z) = 
R v i r (M, z)/c n f w , where R v i T is the cluster's virial radius computed according to the spherical 
collapse model outlined in section 3.1 and appendix A, and c n f w = 5 is the halo concentra- 
tion parameter determined from iV-body simulations. The factor a(M, z) is computed from 
smoothing the (projected) NFW profile using a Gaussian filter of angular smoothing scale 
Og, i-e-, 

a(M, z) = Jo" dx ( x / x g) ex P (-x 2 /x 2 g) f ( x ) ; (52) 

In (1 + c nfw ) - c nfw / (1 + c nfw ) 

where x = 9/9 s and xq = 9q/9 s , with 9 s (M,z) = R S (M , z) / d^z) and d& (z) the angular 
diameter distance to the cluster. The projected NFW profile is encoded in the dimensionless 
surface density profile f(x), which can be found in equation (7) of reference [33]. In our 
analysis we use an angular smoothing scale of 9q = 1 arcmin. 

The mean critical surface mass density S cr (z) is, assuming a flat spatial geometry, given 
by the expression 

V~Xz) = J^ < J dz> dn/dz' X (z) [1-x (z) f x (z>)] , (5.3) 

where x( z ') denotes the comoving radial distance to the redshift z', and (dn/dz')dz' is the 
number density of source galaxies per steradian at redshift (z' , z' + dz'), normalised such that 
iT-bg = Jo dz' dn/dz' gives the source galaxy surface density. As in [3], we assume a galaxy 
redshift distribution of the form 

/ „ \ 2 r / „ \ p\ 

(5.4) 

where for a EuCLiD-like survey we choose f3 = 1, Zq = 0.3, and a source galaxy surface 
density of 7ib g = 40arcmin~ 2 [2]. 

In order for a cluster to be considered detected, its shear signal kg must exceed the 
"noise" of the survey a no \ se by a predetermined amount. Shear detection is limited firstly 
by the intrinsic ellipticity of the background galaxies, and secondly by the number of galaxy 
images lensed by the cluster that fall within the smoothing aperture. Thus, the noise term 
may be estimated as [35] 

anoise = 47r0|n bg ' (5 ' 5) 

where a e denotes the mean dispersion of the galaxy intrinsic ellipticity, and we use a e = 0.1 
in our analysis [2]. 

Defining a signal-to-noise ratio of S/N = 4.5 to be our detection threshold (consistent 
with a survey efficiency of e = 0.75 [36], to be discussed in section 5.3), the expression 

S/N = KG(Mthr( ^ (5.6) 

"noise 

can now be solved for the mass detection threshold M tnr (-z)- This sets a lower limit on the 
cluster mass detectable by lensing at a given redshift z. 
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Figure 1. The left panel shows the division into of the observed number of clusters into 7 redshift 
bins while keeping the cluster count common for all bins. The right panel shows the subsequent 
division of redshift bins i = 1, 7 into 7 mass bins, again with the stipulation that all mass bins contain 
the same number of clusters. 



5.2 Redshift and mass binning 

We consider a survey that observes clusters in the redshift range [z\ ow , Zhigh]- We subdivide 
this range into N z bins in such a way so as to maintain the same number of clusters in all 
bins i in the fiducial cosmology. The resulting bin boundaries 2; m ; n) j and z max ,i then define the 
redshift window functions (4.1). Clusters in each redshift bin are further subdivided according 
to their observed masses into iV m mass bins labelled by j(i), again with the enforcement that 
the number of clusters A^,j(i) should be similar in all bins. 

An immediate consequence of such a binning scheme is a variation of the mass bin 
boundaries M min j^ and M max jU-\ with redshift because of (i) the z-dependence of the mass 
detection threshold M t ^ r (z), and (ii) the rarity of high-mass clusters at high redshifts. For 
the latter point, we impose in practice an absolute high-mass cut-off of Mh; g h = 10 16 M@, i.e., 
the upper limit of the last mass bin, M m3X ^ N u\, is always equal to M^gh a t all redshifts. This 
number also sets the high-redshift cut-off z^igh, which is defined to be the redshift at which 
the mass detection threshold M t ^ r (z) exceeds Mhi g h- The lower cut-off is set at z\ ow = 0.01, 
since the survey contains a negligible number of clusters below this redshift because of the 
small volume and a large detection threshold. 

Figure 1 illustrates the division of the observed number of clusters into redshift and 
mass bins in the case N z = N m = 7. The left panel shows the first division in redshift, while 
the right panel shows the subsequent division of redshift bins i = 1, 7 into mass bins. 

5.3 Completeness and efficiency 

The completeness of a cluster survey is defined as the fraction of clusters actually detected as 
peaks by the cluster finding algorithm, while the efficiency is the fraction of detected peaks 
that correspond to real clusters. In general these quantities can be established precisely only 
with the help of mock cluster catalogues generated from iV-body simulations (see, e.g., [17, 
33, 36-39]). Here, we adopt the same simplistic approach taken in reference [16], and assume 
both c and e to be mass- and redshift-independent. Then the effect of a survey completeness 
and efficiency not equal to unity amounts to increasing the uncertainty on each individual 
data point by a factor -^/[(l/e — 1) + l/e]/c, which can be incorporated into the forecast 
analysis at the level of the likelihood function. We shall return to this point in section 6.2 



where we discuss explicitly the construction of the likelihood function. Suffice it to say for 
now that we adopt the values c = 0.70 and e = 0.75, which may be reasonably expected for 
the LSST [16], and which are likely very conservative when applied to Euclid because of its 
much narrower point spread function. 

5.4 Synthetic data sets 

We summarise here the mock EuCLiD-like data sets we generate and use in our parameter 
sensitivity forecast. 

• A cluster data set in the redshift range z G [0.01, 2hi g h], and the mass range M G 
[Mthr(-z), 10 16 Mq], where M t h r (z) denotes the redshift-dependent mass detection thresh- 
old as described in section 5.1, and Zhi g h is defined as the redshift at which M t -^ r {z) 
exceeds 10 16 M@ as discussed in section 5.2. We slice the redshift- and the mass-space 
into N z and N m bins respectively according to the scheme detailed in section 5.2. 

• Mock data from a PLANCK-like CMB measurement, generated according to the proce- 
dure of [32]. Note that although we do not use real Planck data [40], only synthetic 
CMB data of comparable constraining power, we shall continue to refer to this synthetic 
data set as "Planck data" when discussing parameter constraints. 

• We use also the cosmic shear auto-correlation power spectrum, the galaxy clustering 
auto-spectrum, and the shear-galaxy cross-correlation power spectrum that will be 
derived from a EuCLiD-like photometric survey. The procedure for generating these 
mock data sets has already been described in [3]. We assume an exactly known linear 
galaxy bias for the galaxy clustering power spectrum. 

6 Forecasting 

We now describe our parameter sensitivity forecast for a EuCLiD-like photometric survey 
including a measurement of the cluster mass function. The forecast is based on the con- 
struction of a likelihood function for the mock data, whereby the survey's sensitivities to 
cosmological parameters can be explored using Bayesian inference techniques. 

6.1 Model parameter space 

Reference [3] considered a 7-parameter space spanned by the physical baryon density u>b, the 
physical dark matter (cold dark matter and massive neutrinos) density co dra , the dimension- 
less Hubble parameter h, the amplitude and spectral index of the primordial scalar fluctua- 
tions ^4 S and n s , the reionisation redshift z re , and the neutrino density fraction f u = w v /wd m , 
with u v = ^771^/(94.1 eV). In the present analysis we extend this model parameter space 
to include also the possibility of a non-standard radiation content, quantified by the effec- 
tive number of massless neutrinos N^, as well as three dynamical dark energy parameters 
6^ = (wo, w a ,(%), taking the total number of free parameters to eleven: 

0(H) = (0(8), qQ) = ((oj h ,oj dm ,h,A s ,n s ,z re ,f u ,NZ£),(w ,w a ,c 2 s )) . (6.1) 

As in [3] we assume only one massive neutrino state, so that N^g = 2.046 + AN, where AN 
parameterises any non-standard physics that may induce a non-standard radiation content. 
Note that the 0.046 contribution in N^ comes from non-instantaneous neutrino decoupling 



and finite temperature QED effects, and should in principle be shared between both the 
massless and the massive neutrino states. However, in practice, the precise treatment of this 
small correction has no measurable effect on our parameter forecast. Lastly, we remark that, 
parameterised as such, AN can run from the lowest value of —2.046 to anything positive. 
Many popular models with non-standard radiation contents associate a positive AN with 
additional relativistic particle species such as, e.g., sterile neutrinos [41, 42]. A negative AiV 
can however arise in, e.g., models with extremely low reheating temperatures [43, 44]. 

For the non-dark energy part of the parameter space, our fiducial model is defined by 
the parameter values 

9J§ = (0.022, 0.1126228, 0.7, 2.1 x 10~ 9 , 0.96, 11, 0.00553, 2.046). (6.2) 

For the dark energy sector, we begin with the fiducial values 0g d = (— 1, 0, oo) corresponding 
to dark energy in the form of a cosmological constant. The first part of our analysis (up to 
and including section 7.4) will also be performed with the dark energy sound speed fixed at 
Cg = oo, i.e., homogeneous dark energy. We shall return to dark energy density perturbations 
in section 7.5, and study the constraints on the dark energy sound speed under a variety of 
assumptions for the fiducial dark energy parameter values ©g d . 

6.2 Likelihood function 

Given a theoretical prediction A^h for the observable number of clusters in a specific red- 
shift and mass bin, the probability of actually observing A b s clusters follows a Poisson dis- 
tribution of iVobs degrees of freedom. However, the imperfect completeness and efficiency 
of the survey necessitate that we rescale uncertainty on each data point by an amount 
/ _1 = y/[(l/e — 1) + l/e]/c (see section 5.3). We accomplish this by defining an effective 
number of observed clusters -/V bs = / 2 A b s , and likewise an effective theoretical prediction 
Ath — f 2 ^th- The effective probability distribution is then 



£p( 



N 



N ohs 



A obs |A th ) - : ^ T exp 
AW 



■N t 



th 



(6.3) 



In a real survey, the effective observed number of clusters A b s in any one bin is necessarily 
an integer so that equation (6.3) applies directly. In our forecast, however, A b s corre- 
sponds to the theoretical expectation value of the fiducial model which generally does not 
evaluate to an integer. To circumvent this inconvenience, we generalise the likelihood func- 
tion (6.3) by linearly interpolating the logarithm of the discrete distribution £p in the interval 
[floor (A obs ), ceiling(iVobs)], i.e., 



ln£(A b s |A th ) = 1 + floor(iV b s ) - A obs ln£ P floor(A obs )|A th 

~ n ~ / ( 6 - 4 ) 

+ (AUs - floor (JVobs)) hi£p (ceiling(A obs )|A th ) . 

The total cluster log-likelihood function is then obtained straightforwardly by summing In C 
over all redshift and mass bins. 

7 Results 

7.1 Impact of the number of bins 

We examine first how the parameter sensitivities of the cluster survey depend on the number 
of redshift and mass bins used. Figure 2 shows the posterior standard deviations for a range 
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Figure 2. Dependence of the posterior standard deviations, a, for selected cosmological parameters 
for CMB+clusters on the number of rcdshift bins N z and mass bins N m . All numbers have been 
normalised to the corresponding N z = N m = 1 result. 

of cosmological parameters derived from a combination of synthetic Planck and cluster data 
as functions of N z and N m , normalised to the corresponding N z = N m = 1 result. 

A general trend is immediately clear: for the parameters wq, w a and h, while increasing 
the number of mass bins results in moderate gain, it is the number of redshift bins used that 
contributes mostly to improving the parameter sensitivities. For example, in the case of a 
fixed N z = 1, the number of mass bins needs to be increased to 10 or more in order for the 
sensitivities to improve by 30%, while the same or better sensitivities can be achieved already 
for a fixed N m = 1 simply by splitting the data into three redshift bins. For N*a, ^2 m u 
and co m , the improvements in sensitivity are similar in both directions, and can be traced 
to the fact that these parameters are primarily responsible for the shape and the overall 
normalisation of the cluster mass function, less so the redshift dependence. See section 7.2. 

Some parameter sensitivities continue to improve beyond N z = 10 and one might there- 
fore argue for using a very large number of redshift bins. However, for N z = 10, the narrowest 
bin typically has a width of order Az = 0.03; pushing much beyond N z = 10 may cause 
our results to lose their robustness against redshift uncertainty (see section 4.1). Further- 
more, when cluster data are used in conjunction with angular power spectra from cosmic 
shear and/or galaxy clustering, the gain in going beyond N z = 10 is substantially reduced. 
Henceforth, we shall adopt N z = A^ m = A r b; n = 10. Figure 3 shows the marginalised joint 
two-dimensional posteriors in the (wo,w a )- and the (A^ 1 ,w m )-subspace from CMB+clusters 
for this configuration, as well as for N^ n = 1, 2. 
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Figure 3. Marginalised joint two-dimensional 68% and 95% credible contours from the CMB+clusters 
data set. The default redshift and mass binning configuration for the cluster data is -/Vbin = 10 (light 
red), but we also show the results for 7Vbi n = 1 bin (dark red) and 2 (red). 



Table 1. Posterior standard deviations for the parameters tu m , h, ^m^, N^g, wq, and w a derived 
from various combinations of data sets. Here, "c" denotes Planck CMB data, "g" galaxy auto- 
spectrum (11 redshift bins), "s" shear auto-spectrum (2 bins), "x" shear-galaxy cross-correlation, and 
"cl" the cluster data (10 redshift bins, 10 mass bins). Also shown is the wo~w a figure-of- merit (FoM), 
defined in equation (7.1) as (a(wo)cr(wa))~ 1 ■ 



Data 



1(T x a(uj m ) 100 x a(h) a^m^/eV ct{N^ 1 ) a(w ) a{w a ) FoM 



csgx 


1.36 


ccl 


1.08 


csgxcl 


0.35 



0.858 
0.323 
0.231 



0.022 
0.039 
0.0098 



0.069 
0.031 
0.019 



0.077 
0.038 
0.025 



0.22 59 

0.16 165 
0.085 454 



7.2 Probes of the expansion history versus probes of the power spectrum 

The main advantage of the cluster mass function (with redshift binning) over the shear and 
the galaxy power spectra is that it is highly sensitive to those parameters that govern the 
linear growth function and hence (in the case of standard gravity) the expansion history 
of the universe. This makes redshift-binned cluster measurements especially powerful for 
constraining dark energy parameters, as well as for establishing the reduced matter density 
il m = co m /h 2 . Furthermore, because the normalisation of the cluster abundance is directly 
sensitive to the physical matter density u m , the Hubble parameter h can also be very effec- 
tively constrained. As shown in table 1 and figure 4, CMB+clusters ( "ccl" ) performs up to a 
factor of two better than CMB+shear+galaxies ("csgx") when it comes to placing constraints 
on wo, w a , uj m , and h. 

In contrast, the sum of neutrino masses ^2,rn v is significantly less well measured by 
clusters than by the shear and the galaxy power spectra. This is because firstly, Yl m v pl a Y s 
a negligible role (compared with, e.g., dark energy parameters) in the redshift dependence 
of the late-time linear growth function. Secondly, although the shape of the cluster mass 
function is in principle also subject to a mass-dependent suppression due to neutrino free- 
streaming (e.g., [45, 46]), the actual range of cluster masses probed by a realistic cluster 
survey is very narrow (see figure 1), so that the suppression can be easily be mimicked 
by other effects such as an excess of relativistic energy density or simply a smaller initial 
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Figure 4. Marginalised joint two-dimensional 68% and 95% credible contours from the CMB+clusters 
data set ("ccl", blue), CMB+shear+galaxies ("csgx", green), and all data sets ("csgxcl", black) for 
various parameters, using the default binning configuration of A^in — 10 for the cluster data. 

fluctuation amplitude. 

Interestingly, a non-standard radiation content as parameterised by iV^g , although it has 
no direct effect on the late-time expansion or growth history, is exceptionally well constrained 
by CMB+clusters. This can be understood as follows: using CMB data alone, N^g is strongly 
degenerate with u> m and h. However, because the cluster mass function is directly sensitive 
to u m and h, it very effectively lifts any degeneracy of these parameters with N^ when used 
in combination with CMB data. As shown in the lower right panel of figure 4, very little 
degeneracy remains between JV2* and w m for the CMB+clusters data set. A more telling 
illustration of how the binned cluster data removes the (A r e ff,w m )-degeneracy can be found 
in the right panel of figure 3: Here, when only one redshift and mass bin is used, the cluster 
mass function is primarily sensitive to the fluctuation amplitude on small scales so that the 
(iV e ff,u; m )-degeneracy persists in the CMB+clusters fit. However, as soon as access to the 
linear growth function and some shape information become available through as little as 
N z = N m = 2 bins, the degeneracy becomes broken because of the growth function's direct 
dependence on f2 m and of the normalisation's dependence on u> m . 

7.3 Combining all data sets: constraints on neutrino parameters 

Perhaps the most noteworthy result of table 1 is that, while CMB+shear+galaxies ("csgx") 
and CMB+clusters ("ccl") are well-suited to measuring different parameters and are hence in 
a sense complementary to each other, the combined usage of all data sets, i.e., the "csgxcl" 
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combination, always leads to significant enhancements in all parameter sensitivities. This 
result can be understood from figure 4, where it is clear that the 'csgx" and "ccl" datasets 
give rise to almost orthogonal parameter degeneracy directions. In combination these data 
sets conspire to lift each other's degeneracies. 

For neutrino parameters, it is interesting to note that while "csgx" and "ccl" return 
a(Y^m u ) = 0.022 eV and 0.039 eV respectively, in combination the sensitivity improves 
to a(^2m u ) = 0.0098 eV. This is as good a sensitivity as was found earlier in [3] from 
the "csgx" data set, but for a much simpler 7-parameter cosmological model. Thus, we 
again conclude that Euclid, in combination with Planck CMB data, will be able to probe 
neutrino masses at 5<r precision or better. Likewise, the "csgxcl" data set is now sensitive to 
N£ff at cr(iVig) = 0.019, meaning that for the first time the small deviation of 0.046 from 3 
in the fiducial N£g can be probed with 2<r precision. 

7.4 The dark energy figure-of-merit 

It is useful to quantify the constraining power of an observation (or a combination of obser- 
vations) in terms of a figure-of-merit (FoM). Here, we define a dark energy FoM, 

FoM = (a{w )a(w a ))-\ (7.1) 

which is simply the reciprocal of the product of the posterior standard deviations for wo 
and w a , so that the larger the FoM the better the constraining power. 

To relate our FoM to that adopted in the Euclid Red Book [2] , we note that the latter is 
defined as FoM = (a(w p )a(w a ))~ 1 , following from a parameterisation of the dynamical dark 
energy equation of state of the form w{a) = w p + w a (a p — a), where the "pivot" scale factor a p 
is chosen such that the uncertainties on w p and w a are uncorrelated. Our parameterisation 
(equation (2.1)) corresponds to setting a p = 1. The two parameterisations of w(a) are in 
principle equivalent in the sense that the areas of the resulting error ellipses in the (wq, w a )- 
and (iUp, io )-subspaces are the same [47]. In the case of the (w p , u> a )-parameterisation, be- 
cause of the absence of correlations the product a(w p )a(w a ) represents a reasonable estimate 
of area of the 68%-credible region. In our (wo, w a )-parameterisation, however, strong corre- 
lations exist between wo and w a , as is evident in figure 3, such that the product a{wo)o{w a ) 
always overestimates the area of the error ellipse. It follows that 

{a{w )a{ Wa ))- 1 < {a{w p )a{w a ))-\ (7.2) 

i.e., our FoM defined in equation (7.1) is conservative in comparison with that of [2]. 

The FoMs for various fiducial models and combinations of data sets are presented in 
table 2, where it is immediately clear that the Euclid science objective of reaching a FoM 
in excess of 400 in the context of the ACDM model can be realised using the "csgxcl" data 
combination. We note that this number is formally much smaller than the estimate of 4000 
quoted in [2]. This can be traced to (a) the conservative nature of our FoM definition, (b) the 
larger parameter space used in our analysis (constraints on w p and w a quoted in [2] have 
been obtained assuming fixed values of ^ m v = eV and N£g = 3.046), and (c) the Fisher 
matrix approach adopted in the analysis of [2], which tends to predict constraints that are 
too stringent compared with those from a Monte Carlo analysis — the recent analysis of [48] , 
for example, found a factor of 2 to 3 overestimation of the sensitivity of future cluster surveys 
to wo and w a - 

Lastly, we note that for a given data set combination, there are some variations in the 
sensitivities to wq and w a (especially the latter) depending on our choice of fiducial model. 
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Table 2. Posterior standard deviations for the parameters Wq and w a in various fiducial models 
derived from various data set combinations. See legend in table 1. Also shown is the WQ-w a figure-of- 
merit (FoM). 



Data 


Wo 


W a 


cr(tUo) 


cr(w a ) 


FoM 


ccl 


-1.00 


0.00 


0.038 


0.16 


165 


ccl 


-0.83 


0.00 


0.036 


0.13 


212 


ccl 


-1.17 


0.00 


0.037 


0.20 


133 


ccl 


-1.00 


0.35 


0.030 


0.14 


231 


ccl 


-1.00 


-0.35 


0.041 


0.22 


109 


csgxcl 


-1.00 


0.00 


0.025 


0.085 


454 


csgxcl 


-0.83 


0.00 


0.022 


0.070 


645 


csgxcl 


-1.17 


0.00 


0.024 


0.083 


496 


csgxcl 


-1.00 


0.35 


0.022 


0.070 


658 


csgxcl 


-1.00 


-0.35 


0.025 


0.086 


470 



The FoM likewise varies by up to 40% away from the ACDM value. However, these variations 
across fiducial models for the same data combination are typically minor compared with the 
variations across different data combinations for the same fiducial model. Thus, by and large, 
results for the ACDM model can be treated as representative. 

7.5 Dark energy sound speed and perturbations 

In the analysis so far we have neglected the effect of dark energy perturbations. We now 
introduce dark energy perturbations into the analysis as per the discussion in section 2, and 
investigate the constraining power of a EuCLiD-like survey on this aspect of dark energy. 
We consider four fiducial models differing in their fiducial values of (wo,w a ,c 2 ): 

• Model 1: (w ,w a ) = (-1,0) and c 2 = l 1 (i.e., ACDM) 

• Model 2: (w ,w a ) = (-1.17,0) and cf = 1 

• Model 3: (w ,w a ) = (-0.83,0) and c 2 = 1 

• Model 4: (w ,w a ) = (-0.83,0) and c 2 = 10~ 6 

In all models, we fix w a = because of the pathological behaviour of dark energy pertur- 
bations when crossing the phantom divide w = — 1 [49]. We present only results obtained 
from the "csgxcl" data sets, shown in table 3 and figure 5. All results have been obtained 
assuming a top-hat prior on logc 2 of logc 2 £ [—10,2]. Note that the posteriors for logc 2 
in some models clearly show no peak structure so that logc 2 can only be constrained from 
one side. In such cases, instead of the posterior standard deviation, we quote in table 3 the 
appropriate one-sided limits. 

For the three models with c 2 = 1, the Jeans mass Mj, defined in equation (A. 6), is of 
order 10 23 M at z = and lies well above the maximum observed cluster mass Mhi g h- This 
tells us immediately that these models are undetectable by a EtJCLlD-like cluster survey. 
However, as c 2 drops below ~ 10~ 2 (Mj ~ 10 20 M ), the additional mass-dependence it 
produces on the cluster mass function begins to be visible to Euclid, and we see a sharp 



1 The choice of c^ for the ACDM model is immaterial, since by definition dark energy does not cluster in 
this model and Sq = is automatically implemented in Camb. 



15 



Table 3. Posterior standard deviations for wo and log c^ in various fiducial models derived from the 
combined "csgxcl" data set. See legend in table 1. In those cases where logCg can only be constrained 
from one side, we quote the corresponding 68% (95%) credible limits. 
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Figure 5. Marginalised one-dimensional posterior probability density for logc!? for four different 
fiducial models (see labels in plots) from a "csgxcl" fit. 

decrease in the posterior probability for log cj; in figure 5. In the case of models 2 and 3 (which 
have wq = —0.83, —1.17), the corresponding posterior probabilities eventually drop to zero, 
thus allowing us to place a lower limit on c;?: logc^ > —4.0 (95% CL). No constraints on 
Cg are however available for model 1 (the ACDM model), because dark energy perturbations 
are generally suppressed by a factor of 1 + wo relative to dark matter perturbations and are 
hence practically nonexistent in the vicinity of wq = —1. 

The remaining model 4 with Cg = 10~ 6 has a remarkably well-constrained dark energy 
sound speed, a(logCg) = 0.13. This is a consequence of Mj (~ 4 x 10 14 M & at z = 0) lying 
within the cluster mass range probed by Euclid, so that the transition from minimal dark 
energy clustering to full dark energy clustering can be very effectively observed with the 
help of cluster mass binning. We emphasise however that our choice of the fiducial value of 
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Wo = —0.83 also plays a crucial role towards establishing this remarkable sensitivity, again 
because of the 1 + u>o suppression suffered by the dark energy perturbations relative to the 
dark matter perturbations. 

8 Conclusions 

In this paper we have considered the constraining power of a EuCLiD-like galaxy survey on 
cosmological parameters in conjunction with Planck CMB data. This study is an extension 
of our previous investigation in [3] , in that we have included in the present analysis mock data 
from the Euclid cluster survey in addition to the angular cosmic shear and galaxy power 
spectra expected from the photometric redshift survey, and we have expanded the parameter 
space to encompass also dynamical dark energy as well as the possibility of (small-scale) dark 
energy perturbations. 

We find that, comparing CMB+clusters and CMB+shear+galaxies, the former has up 
to a factor of two better sensitivities for those parameters that affect only the late-time growth 
and expansion history of the universe, i.e., those parameters that determine the dynamical 
dark energy equation of state, and the Hubble parameter. This improvement arises chiefly 
from our adoption of redshift binning for the observed clusters, which allows us in particular 
to probe the transition from matter to dark energy domination. Neutrino masses, on the 
other hand, are not particularly well-constrained by CMB+clusters (cr(^2 m u ) = 0.039 eV in 
a 10-parameter model), clearly because they do not play a major role in the overall linear 
growth of matter density perturbations. Importantly, however, the degeneracy directions 
of CMB+clusters and CMB+shear+galaxies are largely orthogonal. This means that even 
though neither data set performs particularly impressively for any one cosmological param- 
eter, in combination they help to lift each other's degeneracies. The sensitivities to Y2 rn v 
from CMB+shear+galaxies+clusters is, for example, cr(Yl m u ) = 0.0098 eV in a 10-parameter 
model, which is as good as that obtained previously in [3] from CMB+shear+galaxies for a 
much simpler 7-parameter ACDM model. Thus, we can conclude again that a EuCLiD-like 
survey has the potential to measure neutrino masses at 5<r precision or more. 

For the dark energy parameters, we find that the combination of CMB+shear+galaxies+ 
clusters results in a dark energy figure-of-merit (FoM), defined in this work as (a(wo)a(w a ))~ 1 , 
of 454 for a ACDM fiducial cosmology, with upward variations of up to 45% for fiducial cos- 
mologies in which wo ^ — 1 and w a ^ 0. Although this number is clearly sufficient to meet 
the Euclid science objective of reaching an FoM beyond 400, we note that it is significantly 
smaller than the estimate of 4000 quoted in the Euclid Red Book [2] . This discrepancy can 
be traced to our more conservative definition of the FoM, as well as the smaller parameter 
space and the Fisher matrix forecast method (which tends to be too optimistic compared 
with a Monte Carlo forecast) adopted in the analysis of [2]. 

Finally, we investigate the detectability of dark energy perturbations, parameterised in 
terms of a non-adiabatic fluid sound speed Cg . Along the way we also introduce a model of the 
cluster mass function that incorporates the effects of Cg based on solving and interpolating 
the spherical top-hat collapse in the known limits of Cg — > oo (homogeneous dark energy) and 
Cg = (dark energy comoving with nonrelativistic matter). We find that for values of the 
dark energy sound speed whereby the associated Jeans mass lies within the mass detection 
range of the cluster survey, dark energy perturbations imprint a distinct step-like signature 
in the observed cluster mass function. With the help of cluster mass binning, this signature 
makes these models distinguishable from those in which the Jeans mass lies well outside (both 
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below and above) the detection range. Indeed, using 10 mass bins per redshift bin, we find 
that logCg can be measured to a la precision of 0.13 in these models, as long as the fiducial 
value of wo deviates from —1 by as much as is presently allowed by observations. 

Note added: As this work was in its final stage of completion, we learnt of the in- 
vestigation of [50] which considered the Euclid cluster survey's sensitivities to neutrino 
parameters. While a full comparison is difficult because of the generally different assump- 
tions about the survey parameters and the model parameter space, where the assumptions 
do to some extent coincide the two analyses appear to be compatible. 
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A Interpolating the spherical collapse between the two limits of c 



The spherical top-hat collapse model is exactly defined only in the limits c 2 — > oo and Cg = 0. 
In the first case, the dark energy component is non-clustering. The dark matter and baryon 
components alone suffer gravitational collapse, so that the overdense region preserves its 
top-hat density profile throughout the collapse, with a comoving radius X given by 

X+ n X= ~^Ga 2 p m 5 m , (A.l) 



Sm(r) = [l + 6 m (Ti)] 



1 (A.2) 



where 

~ X(n) ~ 

Mr). 

follows from conservation of the total mass of nonrelativistic matter M m in the top-hat 
region. The virial radius R v i T , defined as the physical radius of the top-hat at the moment 
the collapsing region fulfils the virial theorem, is as usual one half of physical radius at turn- 
around, and the virial mass M v - ir is identical to M m , which we also equate with the mass of 
the cluster M. 

In the second case, an exactly vanishing non-adiabatic dark energy sound speed means 
that, like nonrelativistic matter, dark energy density perturbations also evolve identically on 
all scales. This again leads to the preservation of the top-hat density profile and hence the 
conservation of M m in the region defined by the comoving radius X, now determined by 

§ + H ^ = -4vrGa 2 [p m 5 m + p Q 8 Q ] . ( A.3) 

A conservation law can likewise be written down for the clustered dark energy component 

- 0, (A.4) 

where pq denotes the dark energy density in the top-hat region. We assume that at any one 
time the clustered dark energy contributes a mass [28] 

M Q (r) = ^p Q ( T )5 Q (T)[a(r)X(r)f (A.5) 
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to the total mass of the system. This clustered dark energy takes part in virialisation, defined 
here as the instant the system satisfies the condition d 2 I/dt 2 = 0, where I = (2/5) (M m + 
MQ)(aX) 2 is the top-hat's moment of inertia, and t is the cosmic time (dt = adr). The 



physical radius of the top-hat at this instant is the virial radius R v 



aX[T v 



and the 



total mass the virial mass M v ; r = M m +Mg(r v i r ). Linearised forms of equations (A.l) to (A. 4) 
can be found in references [9, 28]. 

Between the two limits, any finite nonvanishing c 2 necessarily causes the overdense 
region to evolve away from the top-hat configuration and to become ill-defined. The absence 
of strict conservation laws in this intermediate regime also renders the system not readily 
soluble. Nonetheless, the transition between the two known limits have been investigated 
using a quasi-linear approach in [9, 29], where it was found that, for a fixed c 2 and at 
a given collapse redshift z, the transition generically results in a step- like feature in the 
linear threshold density <5 C , as well as in the quantities A v j r = 3M m /47rp m (r v i r )i?^ ir and 
r/ v j r = Mg(r v i r )/M m , as M m is varied from M m <C Mj to M m 3> Mj, with 



Mj(a) 



4tt 



Pra{a) 



oAj(a) 



9.7 x 10 23 c 3 Q m a' 3 (n m a^ 3 + n e 



-3f« da iil +W (ai)]/aiy 3 / 2 h -l MQ (A6) 



denoting the Jeans mass, and Aj = 2irc s /T-L the corresponding comoving Jeans length. 

To incorporate this feature in the cluster mass function for our MCMC analysis, our 
strategy is as follows. 

1. For each given set of cosmological parameters, we solve the spherical collapse model for 
S c , A v i r , and ?7 v i r in the two limits of c 2 using equations (A.l) to (A. 4) as functions of 
the collapse redshift z. 

2. At each redshift, we interpolate the two limits using the formulae 



S c (M m ) 
A vir (M n 

»7vir(M m 



A(<£°,<£)tanh 



A(A^ r ,A° ir )tanh 



Ay log M m - log 



Mj 



A 2 I log M m - log 



oo X Q\ 



+ S(C«5, 

Mj 



B 2 



(A.7) 
+ S(A-,A° ir ),(A.8) 



A(rf. 



~,^ ir )tanh 



A 3 logM n 



i A M 
" l0g W 



+ ^'( 7 ?vir' r /vir)' 



(A. 



where x°° denotes the value of x = S c , A v i r , ?y v i r in the Cg — > oo limit, x u the c 2 = limit, 

-°°,2; ) = [x 



50 + x°)/2 represent their difference 



and A(x°°,x°) = {x°° - x°)/2 and S(x ( 
and sum respectively. The parameters -Ai,2,3 an d ^1,2,3 are fitting coefficients, adjusted 
to fit respectively the sharpness and location of the transition. At z = and for 
(wo,w a ) = (—0.8,0), we find that ^1,2,3 — 1-4, and .#1,2,3 — 20 reproduce the quasi- 
linear results of [9, 29] in the region immediately below the Jeans mass quite well. For 
simplicity, we adopt these parameter values for all redshifts and cosmological models. 



3. Lastly, we identify the virial mass with the cluster mass, i.e., M v i r = 
the functions 5 C (M, z) and R v i T (M, z) according to 

5 C (M, z) = 5 c (M vir (M m ), z) = 6 c (M m , z), 
R vir (M, z) = R vh (M vh (M m ), z) = R vir (M m , z), 



M, and construct 

(A.10) 
(A.ll) 
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which are used in sections 3.1 and 5.1 determine the cluster mass function and the mass 
detection threshold, respectively. 
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